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Abstract 

Using tools from spectral analysis, singular and regular perturbation theory, we develop a systematic 
method for analytically computing the approximate price of a large class derivative-assets. The payoff of 
the derivative-assets may be path-dependent. Additionally, the process underlying the derivatives may 
exhibit killing (i.e., jump to default) as well as combined local/nonlocal stochastic volatility. The nonlocal 
component of volatility may be multiscale, in the sense that it may be driven by one fast-varying and 
one slow-varying factor. The flexibility of our modeling framework is contrasted by the simplicity of our 
method. We reduce the derivative pricing problem to that of solving a single eigenvalue equation. Once 
the eigenvalue equation is solved, the approximate price of a derivative can be calculated formulaically. 
To illustrate our method, we calculate the approximate price of three derivative-assets: a vanilla option 
on a defaultable stock, a path-dependent option on a non-defaultable stock, and a bond in a short-rate 
model. 
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1 Introduction 



The spectral represent ation for the tr ansition density of a general one-dimensional diffusion was obtained 
in a seminal paper by 



McKeanl (| 1956!) . Since that time, spectral theory - and more specifically, the study 
of eigenfunction expansions of linear operators - has become an essential tool for analysing diffusions. As a 
diffusion often serves as the underlying process on which financial models are built, it is not surprising that 
methods from spectral theory have made their way into mathematical finance as well. 

In particular, many problems related to the pricing of derivative-assets have been solved analytically by 
using methods from spectral theory. An overview of the spectral method applied to derivative pricing is 
as follows. Using risk-neutral pricing, one expresses the value of a derivative-asset u(t, x) as a risk-neutral 
expectation of some function of the future value of an underlying process X . Mathematically, this is expressed 
as 

(1.1) u(t,x)= : E x [H(X t )] = J H(y)p(t,x,y)dy. 

Here, p(t, x, y) is the transition density of the X under P. If it turns out that the ininitesmal generator 
£ of the underlying process is self-adjoint on a Hilbert space with weighting measure m(x)dx and if the 
spectrum of L is purely discrete, then the transition density of X has an eigenfunction expansion 

(1.2) p(t 7 x,y) = m{v)^e- Xnt ip n {y)ip n {x), 

a 

where {A n } are the eigenvalues of (— L) and {ipn} are the corresponding eigenfunctions 

— L ip n = X n ip n . 

The value of a derivative-asset can then be expressed analytically by inserting (|1.2|) into (|l.ip 
u(t,x) = ^2c n e~ Xnt ip n (x), c n = (tp n ,H) := / H(y)ip n (y)m(y)dy. 

Under some basic assumptions, the infinitesimal generator of a general one-dimensional diffusion 

(1.3) L = ]^a 2 {x)d 2 xx + b(x)d x - k(x), x £ (e 1 ,e 2 ), 

with domain dom(£) (described in appendix IA.2[) is always self-adjoint on the Hilbert space 'M = L 2 (I,m), 

where / C K is an interval with endpoints e\ and ei and m is the speed density of the diffusion 

2 ^ / r 2b(y) i 
a 2 {x) \J X0 a 2 (y) 



(1.4) m(x) := — - exp ( / - dy ) . (speed density) 



1 An operator L is self-adjoint on a Hilbert space Jf with inner product (-, ■) if dom(£) = dom(£*) and (£/, g) = (/, £g) for 
au />9 € dom(£). Please see appendix I A. 1 1 for a brief review of self-adjoint operators in Hilbert Spaces. 
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The lower limit of integration xq £ I is arbitrary. Thus, when a one-dimensional diffusion is adequate for 
describing the dynamics of an underlying, the spectral method outlined above serves as a powerful tool for 
analytically pricing derivatives on that underlying. Among the topics that have been addressed by applying 
spectral methods to one-dimensional diffusions are option pricing (b oth vanilla and exotic), mortgages valu 



ation, i nterest rate modeling, volatility modeling, and c r edit risk (see 



(2002); 



2004); 



Davvdov and Linetskvi (1 2003) ; 



Linetskvi (l2004bl) ; 



Albanesc, Campoli eti. Carr. and Lipton (12001); 



Davvdov and Linetskv (2001) 



Albanese and Lawi (200 



Linetskv 



Albanese and Kuznetsov 



Lewis! dl998 HLipton and McGheeJ (12 002) ; 



Goldstein and Keirstead 



,1997); 
2006)). 



Gorovoi and Linetskv 



2004); 



Gorovov and Linetskv 



20071) : 



Carr and Linetskv 



2006) 



Linetskv 



2004a. c 



A use ful reference on the topic of spectral methods for one-dimensional diffusions in finance is 



Linetskvi (|2007f ) 



As widely applicable as one-dimensional diffusions are in finance, there are applications in which one- 
dimensional diffusions are not adequate for describing the dynamics of an underlying. This is the case, for 
example, in a stochastic volatility setting, where the volatility of the asset that underlies a derivative is con- 
trolled by (possibly multiple) nonlocal diffusions. Ideally, one would like to employ techniques from spectral 
theory to solve problems that relate to multidimensional diffusions. Unfortunately, whereas the infinitesimal 
generator of a one-dimensional diffusion is practically guaranteed to be self-adjoint, the infinitesimal gen- 
erator of a multidimensional diffusion is only self-adjoint when the drift vector satisfies certain constraints 
imposed by the vol atility matrix. The drift constrai n t is no t satisfied by any of the most prominent stochastic 



volatility models 



Hestonl ( 199 
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Hull and White! (|l987l) . 



Hagan. Kumar. Lesniewski. and Woodwax 



Recently, 



dl ( 2002^ 



Stein and Steinl f|199lh and the SABR model by 



which complicates the use of spectral methods. 



Fouque. Jaimungal. and Lorigl ([201 if ), show one way to deal with this issue. By combining tech- 



niques from singular perturbation theory and spectral theory, the authors are able to express the approximate 
price of a (possibly path-dependent) option as an eigenfunction expansion, even though the infinitesimal gen- 
erator of the two-dimensional diffusion t h ey wo rk with is not self-adjoint. As notable as their work is, the 



results of 



Fouque. Jaimungal. and Lorigl (|2011l) are valid only when the asset underlying the option is a 



Black-Scholes-like geometric Brownian motion (GBM) with fast mean-reverti ng stochastic volatility. 



In this paper, we extend the work of 



Fouque. Jaimungal. and Lorigl ([201 11 ) in four important ways. 



1. As a "base" model, we work with a general one-d i mensi onal diffusion dX t = v(X t )dt + a(X t )dWt 



This is in contrast to 



Fouque. Jaimungal. and Lorigl f|201 lh , where the only base model considered is a 



GBM: dX t = fiX t dt + aX t dW t . 
2. The general diffusion we work with may exhibit killing (jump to default) at a rate h{X t ) > 0. In the 
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GBM case considered in lFouque. Jaimungal. and Lorigj ([201 II ). X is always strictly positive 



3. To our general diffusion we add two factors of nonlocal volatility: a(X t ) — > a(X t )f(Y t , Z t ). The first 
factor Y is a fast- varying factor. The second factor Z is slow- varying. Thus, our model i s a m ultiscale 



stochastic volatility model. Again, this is in contrast to 



Fouaue. Jaimungal. and Lorigl (|201ll) . where 



the analysis is limited to a single fast mean-reverting factor of volatility aX t — > f(Y t )X t . 

4. In changing from the physical probability measure to the risk-neutral pricing measure, we consider a 
class of market prices of ri sk that is general enough to treat cred it, equity, and interest rate derivatives 



in a single framework. In 



Fouaue. Jaimungal. and Lorigl (| 2 llh the form chosen for the market price 



of risk restricts the authors to equity derivatives only. 



As 



Fouaue. Jaimungal. and Lorigl (|201l[ ). we will derive an eigenfunction expansion for the approximate 



price of a derivativ e-asset despite the fact that th e infin itesimal generator we consider is not (in general) 



self-adjoint. Unlike 



Fouaue. Jaimungal. and Lorigl ([201 11 ). because our multidimensional diffusion contains 



both a fast- varying and a slow- varying factor of volatili ty, we must combine techniques from both singular 



Fouaue. Jaimungal. and Lorigl (|201 lh , only singular 



and regular perturbation theory to achieve our result. In . 
perturbation techniques are required, due to the presence of a single fast mean-reverting factor of volatility. 

Of course, the idea of combining singular and regular perturbation techniques in a multiscale stochas- 
tic volatility setting is not particularly new or unique. The seminal paper on the sub ject, applied in a 



Black-Scholes-like GBM setting, is due to 



Fouaue. Papanicolaou. Sircar, and Spina 



tionof the singular and regular perturbation methods developed in 



( 20041) . Further applica- 



Fouaue. Papanicolaou. Sircar, and Spina 



(|2004l ) led to papers conc erning bond-pricing, interest rate deriv a tives, credit derivatives, and option pricing 



i n a CEV-like setting (see 



(|200J); 



DeSantiago. Fouaue. and Spina (2008); 



Fouque. Sircar. 



and Spina 



There is also a book by 



2006); 



Cotton. Fouaue. Papanicolaou, and Sircar 



Fouaue, Wignall. and Zho 



Fouaue. Papanicolaou. Sircar, and Solna 



u UoOsL 



Choi. Fouque. and Kim 



201oL 



201ll ). which contains the many of the key 



results from the above mentioned publications. What this paper contributes to the existing literature on mul- 
tiscale diffusions is flexibility and simplicity. From a flexibility standpoint, the meth ods developed in this pa- 



per ar e able to encapsulate, in a unified framework 




Cotton. Fouaue. Papanicolaou, and Sircarl (12 004): 



man y of the results contained in 



Choi. Fouq 



DeSantiago. Fouaue. and Spinal (12 008) ; 



ue. and Kim 



Fouque. Sircar, and Spina 



Fouque. Papanicolaou. Sircar, and Solna 



20041 ): 



Fouque. Papanicolaou. Sircar, and Solna 



2QllD; 



Fouque. Wignall. ai 



20081 ). as well as further results, which are not contained in these works (e.g., jump to default CEV with 



multiscale stochastic volatility, see section 143)) . With regards to simplicity, the spectral method we develop 
reduces the derivative pricing problem to that of solving a single, one-dimensional eigenvalue equation. Once 
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this equation is solved, the approximate price of a derivative-asset can be calculate d formulaically by com- 



puting a few simple inner products. This is in con t rast t o the methods developed in 




Cotton. Fouaue. Papanicolaou, and Sircarl (12 004) ; 



Choi. Fouque. and Kim 



DcSantiag o. Fouaue. and Spinal (12 008') ; 



Fouaue. Sircar, and Sglna 



Fouaue. Papanicolaou. Sircar, and S0lna 



2004|); 



Fouaue. Papanicolaou. Sircar, and Solna 



20111) : 



Fouaue. Wignall. ai 



2008), where, in order to express the approximate price of a derivative- asset, an inhomogeneous partial dif- 



ferential equation (PDE) must be solved. 

The rest of this paper proceeds as follows. In section [5] we introduce a class of models described by 
multiscale diffusions. We also explain the kind of derivative-asset we wish to consider. In section [3] we 
solve (approximately), the problem of pricing a derivative- asset. This is done in several steps. First, using 
risk-neutral pricing, we derive a Cauchy problem which, if solved, would yield the exact value of a derivative- 
asset. Next, we use techniques from singular and regular perturbation theory to formally derive three simpler 
Cauchy problems, which, if solved, would yield the approximate value of a derivative- asset. Finally, using 
eigenfunction expansion techniques, we solve these Cauchy problems explicitly. The solutions are given in 
Theorems 13.11 13.21 and 13.31 In section SI we illustrate our method of pricing derivative-assets with three 
examples. We also provide an appendix, which contains some mathematical results that we use throughout 
this paper. 



2 A Class of Multiscale Models 

Let (CI, 5F, P) be a probability space supporting correlated Brownian motions (W x , W y , W z ) and an exponen- 
tial random variable £ ~ Exp(l), which is independent of (W x ,W V ,W Z ). We shall consider a three-factor 
economy described by a time- homogenous, continuous-time Markov process X = (X, Y, Z), which takes val- 
ues in some state space E = I x R x M. Here, I is an interval in R with cndpoints e\ and e 2 such that 
— oo < ei < e 2 < oo. We assume that X starts in E and is instantaneously killed (sent to an isolated 
cemetery state A) as soon as X leaves /. Specifically, the dynamics of X under the physical measure P are 
as follows: 



(X t ,Y u Z t ) Tl >t 

T I =mi{t>0:Xt^I} 1 

A Tl < t 
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where (X, Y, Z) are given by 



dX t 



dY t 



dZ t 



v(X t )dt + a{X t )f{Y tl Z t )dW t 
-a(Y t )dt+^p(Y t )dW t y , 

6c(Z t )dt + VSg{Z t )dW t z , 



t 1 



(2.1) 



< 



d(W x ,W y ) t 



d(W x ,W z ) t 



Pxz dt, 



d(W y ,W z ) t 



Pyz dt, 



(X , Yq, z ) 



(x,y,z) £ E. 



Here, (p xy , p xz , p yz ) satisfy \p xy \, \Pxz\, \pyz\ < 1 and 1 + 2p xy p xz p yz - p 2 xy - p 2 xz - p* z > so that the 
correlation matrix of the Brownian motions is positive-semidefinite. 

The process X could represent a variety of things. For example, it could represent the price of a stock, 
the value of an index, the risk-free short-rate of interest, etc. More generally, X could represent an exogenous 
factor that controls the value of any or all of the items mentioned above. Under the physical measure P, the 
process X has instantaneous drift v{Xt) and stochastic volatility a(Xt)f(Yt, Zt) > 0, which contains both a 
local component a(X t ) and nonlocal component f(Y t , Z t ). The nonlocal component of volatility f(Y t , Z t ) is 
controlled by two factors: Y and Z. We note that the infinitesimal generators of Y and Z 



are scaled by factors 1/e and S respectively. Thus, Y and Z have intrinsic time-scales e > and 1/5 > 0. 
We assume e << 1 and S << 1 so that the intrinsic time-scale of Y is small and the intrinsic time-scale of Z 
is large. Hence, Y represents a fast- varying factor of volatility and Z represents a slow- varying factor. Note 
that Ly and L s z have the form (| 1 . 3[) with k(x) — for all x £ I. Throughout this paper, we will assume 
that the domain of any operator of the form (|1.3[) is given by equation (|A.3[) of appendix IA.2I 

We are interested in pricing a (possibly defaultable) derivative-asset, whose payoff at time t > may 
depend on the path of X. Specifically, we shall consider payoffs of the form 



Here, r is a random time, which represents the default time of the derivative-asset. Because we are interested 
in pricing derivatives, we must specify the dynamics of {X, Y, Z) under the risk-neutral pricing measure, which 



(2.2) 




(2.3) 



Payoff = H{X t )l {T>t} . 



G 



ia(F t ) - ^=f3(Y t )A(Y t , Z t )j dt + ^PS) dWf , 



we denote as P. We have the following risk-neutral dynamics 

dX t = (b(X t ) - a(X t )f(Y t , Z t )Q(Y t , Z, > ) ,// + >, ( .V, ) / i Y . Z, ) < /I V/ . 
dY t 

dZ t = (5c{Z t ) - V&g{Z t )T{Y u Z t )^ dt + V5g{Z t ) dW t z , 

d(w x ,wy) t = Pxy dt, 

Pxz dt, 

Pyz dt, 

(X ,Y ,Z ) = (x,y,z) € E, 



(2.4) 



d(w x ,w z i 



d(W v ,W z 



where 



dW t x := dW t x 



v{X t ) - b(X t ) 



+ fl(Y u Z t ) dt, 



a{Xt)f(Y t ,Z t ) 
dWf := dW% + A(Y t ,Z t ) dt, 
dW t z := dW t z + T(Y t ,Z t ) dt, 

are driftless BM's under P. We assume (|2.4j) has a unique strong solution. 

As mentioned above, the random time r represents the default time of the derivative-asset. In our 
framework, default can occur in one of two ways. Either default occurs when X exits the interval /, or default 
occurs at a random time t>,, which is controlled by an instantaneous hazard rate h(X t ) > 0. Mathematically, 
we express the default time r as follows 



(2.5) 



T = Tj A T h , 

ti = inf {t > : X t <£ 1} , 



r h = inf ji > : h(X s ) ds > z\ , £ - Exp(l), £ _U_ (X, Y, Z). 



Note that the exponentially distributed r a ndom variable £ is independent of (X, Y, Z) 



Following 



Elliott. Jeanblanc. and Yorl (|2000l ). to keep track of Th, we introduce the indicator process 



Dt '■= I{t>r h }- Denote by D = {D t , t > 0} the filtration generated by D and by F = {£F t) t > 0} the filtration 
generated by {W x , W y , W z ). Define the enlarged filtration G = {5 t ,t > 0} where % t = £F 4 V T) t . Note that 
(X, Y, Z) is adapted to G and r is a G-stopping time (i.e., {r < t} £ St for every t > 0). 

We shall assume our economy includes a risk-free asset, which grows instantaneously at short-rate r(X t ) > 
0. Thus, if our economy includes, for example, a non-dividend-paying defaultable asset S, whose price process 
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is described by St = I{ T>t }X t , where the state space of X was I — (0,oo), then the discounted asset price 
{e~ fo r ( x *) ds S t , t > 0} must be a (P, G)-martmgale. The martingale property can be achieved by setting 
b(X t ) = [r{X t ) +h(X t )]X t and Q(Y t ,Z t ) = in dHJ). The reason for adding the hazard rate h(X t ) to 
the risk-free rate of inter est r(X t ) in the drift of X is to compensate for the possibility of a default (see 



Carr and Linetskvl (j2006f ). Section 2). 



On the other hand, if X only describes the risk-free rate of interest through r(X t ), then in changing from 
the physical measure P to the pricing measure P, one may not have a reason to change the drift of X from 
v{Xt) to b(X t ). However, one may still wish to consider the effect of including a market price of risk. In 
this case, one could set b(X t ) = v{X t ) and keep fl(Y t , Z t ) ^ in (|2.4p . 

We have now described our economy under both the physical and risk-neutral pricing measures, and we 
have specified the kind of derivative-asset we wish to price. However, we have not been specific about certain 
technical assumptions, which we shall need in order to prove the accuracy of our pricing approximation. 
Specific model assumptions can be found in Appendix IA.3I 



3 Derivative Pricing 



We wish to price a derivative- asset whose payoff is of the form (|2.3[) , where the default time r is given by (|2.5[) . 
Using risk-neutral pricing and the Markov property of X, the value u e ' S (t,x,y 7 z) of such a derivative-asset 
at time zero is given by 

rt 



u c ' (t,x,y,z) = E XtVy 



exp 



r(X s )ds)H(X t )l {t>T} 



where (x,y,z) £ E represents th e starting point of the process (X,Y,Z). By conditioning on the path of 



X (see p. 225 of 



Linetskvl (|2007l )) and by using the Feynman-Kac formula, one can show that u e ' S (t, x, y, z) 



satisfies the following Cauchy problem 



(3.1) 
(3.2) 



u t - 5 (0,x,y,z)^H{x) 1 



(x,y,z) G E,te M+, 
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where the operator L e ' 5 is given by 



(3.3) 



C e ' s = -L a + -^=L X +L 2 + \ -M 3 + V^Mi + SM 



L = -p 2 (y)d 2 yy + a(y)d y , 



(3.4) Li = /3(y)(p xy a(x)f(y,z)d x - A(y,zf)d y , 

L 2 = \a 2 (x)f 2 {y,z)d 2 xx + (b(x) - a(x)n(y, z)f(y, z))d x - k(x), 
M 3 = Pxz (3(y)g(z)d 2 z , 
M i = 9{ z ){pxza(x)f(y,z)d x - r(y, z))c\, 
M 2 = ±g 2 (z)d 2 zz +c(z)d z , 
k(x) = r(x) + h(x). 

Aside from the initial condition (|3.2I) . the function u £,l5 (t, x, y, z) must satisfy additional boundary conditions 
(BCs) at the endpoints e\ and e 2 of the interval /. The BCs at e\ and e 2 are understood to be contained 
in the domain of L e ' S and will depend on the nature of the process X near the endpoints of /. Appropriate 
BCs are discussed in appendix IA.21 

From equation (|2.2[) we see that £o = £y- We assume that a diffusion with generator Ly nas an 
invariant distribution n with density tt. In section 13. 1[ it will be important to note that the operator L$ 
with dom(£o) = L 2 (M.,tt) is self-adjoint acting on the Hilbert space L 2 (R,7r). 



3.1 Formal Asymptotic Analysis 

We wish to solve Cauchy problem (|3.1l) - (|3.2p . For general (f, a, f3, A,c, g,T), no analytic solution exists. 
However, we notice that, for fixed 5, the terms in (|3.3p containing e are diverging in the small-e limit, giving 
rise to a singular perturbation. Meanwhile, for fixed e, the terms containing S are small in the small-5 
limit, giving rise to a regular perturbation. Thus, the small-e and small-<5 regime gives rise to a combined 
singular- regular perturbation about the 0(1) operator L 2 . This suggests that we seek an asymptotic solution 
to Cauchy problem (|3.1j) - (|3.2[) . To this end, we expand u e ' S in powers of y/e and VS as follows 

j>0 i>0 

Our goal will be to find an approximation of the price u £ ' <5 « uo.o + \/ewi,o + \/Suo,i- The choice of expanding 
in half-integer powers of e and 5 is natural given the form of L e,s . We will justify this expansion when we 
prove the accuracy of our pricing approximation in Theorem [331 
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Because we are performing a dual expansion in half-integer powers of e and S, we must decide which 
of these parameters we will expand in first. We choose to perform a regular perturbation expansion with 
respect to S first. Then, within each of the equations that result from the regular perturbation analysis, we 
will perform a singular perturbation expansion with respect to e. [^| 



Regular Perturbation Analysis of Equation (|3 . 1 1) 

The regular perturbation expansion proceeds by separating terms in L e ' S and u e ' S by powers of \5 
(3.5) L € ' s =L e + V5M e + SM 2 , = ^ (VPj 3 u), 



j>0 



where 



(3.6) L e = -L + -pLi +L 2 , M e =-=M 3 + M lj 



(3.7) u) = £ (V~e) 



i>0 

Inserting expansions (|3.5[) into PDE (|3.1I) and collecting terms of like-powers of ^/S we find that the lowest 
order equations of the regular perturbation expansion are 

(3.8) 0(1): 0=(-d t + £ e K, 

(3.9) 0(Vd): = (~d t + L e )u{+M t u e Q . 

Now, within equations (|3.8I) and (|3.9p . we will perform a singular perturbation expansion with respect to 
the parameter e. We begin with (I3.8[) . the 0(1) equation. 

Singular Perturbation Analysis of Equation (|3.8[) 

We insert expansions (|3.6p and (I3.7[) into (I3.8[) and collect terms of like-powers of y/e. The resulting order 
(1/e) and (1/y/e) equations are 

(3.10) 0(l/e): = L u , , 

(3.11) 0(l/Ve): O = £o«i,o+-Ci«o,o. 

We note that all terms in Lq and Z\ take derivatives with respect to y. Therefore, if 1*0,0 and ui.o are 
independent of y, equations (13. 10)) and (13. lip will be satisfied. Thus, we choose uq,o = Uo t o(t,x, z) and 



2 Note that we do not take a limit as e and S go to zero simultaneously. 
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ui.q — u\_ o(t, x, z). Continuing the asymptotic analysis, the order (1) and (y/e) equations are 

(3.12) 0(1): = £ u 2 ,Q + (-d t + £ 2 )u ,o, 

(3.13) 0(Ve) : = L u 3fi + £1^2,0 + (-d t + £ 2 ) tti.o- 

where we have used £11/1,0 = in (I3.12[) . Equations (13. 12)) and (|3.13[) are Poisson equations of the form 

(3.14) = L u + X - 

Recall that £0 is a self-adjoint operator acting on L 2 (M,7r). By the Fredholm alternative Q in order for 
equations of the form (|3.14p to admit solutions u G dom(£o) — L 2 (R,ir), the following centering condition 
condition must be satisfied 

(3-15) ( x ) ■= J x (y)<y)dy = o, 

where we have introduced the notation (•) to indicate averaging over the invariant distribution n. In equations 
(|3.12p and (|3.13[) centering condition (|3.15|) corresponds to 

(3.16) 0(1): = (-d t + {L 2 ))u 0i0 , 

(3.17) O(v^): = <£iu 2l o> + (-ft + (-Ca))«i,o. 

The operator (£ 2 ) is given by 

(3.18) (L 2 ) = ^f 2 a 2 (x)d* x +(b(x)-Jla(x))d x -k(x), x € (ei,e 3 ), 
where we have defined 

a 2 (z):=(f(;z)), Jl(z) :=</(•, z)Q(;z)), 

We assume ct 2 (z) < 00 and fii(z) < 00. Given appropriate BCs at e% and e 2 , one can find a unique solution 
Uo,o to PDE (|3.16[) . However, in order to make use of p.!7[) we need an expression for (£iU 2> o). To this end, 
we note from (|3. 121) that 

£qu 2 .o = - (-d t + £ 2 ) w ,o 

= - {-dt + £2) u fi + (-d t + (£ 2 )) «o,o 

\a 2 (/ 2 - a 2 ) 8 2 x - o(/fi - JU) d x ] uo,o. 



J Please refer to Appendix IA. 41 for an discussion of the Fredholm alternative 
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Now, we introduce (j)(y, z) and rj{y, z) as the solutions to the following Poisson equations 

(3.19) £ o = / 2 -a 2 , £ r) = fQ-Jn. 
Using (|3.19p . we can express u 2fi as 

(3.20) w 2 , - - Qa 2 <f> d 2 xx - a r, u ,o + C. 

Note that C is a constant that is independent of y. Now, inserting (|3.4[) and (|3.20p into (£o w 2,o) we find 

(3.21) (<Ci«2,o) = - ((p(pxyafd x - Ajdy^j (^a 2 4>d 2 xx - a-qd^j = -Au Qfi . 
The operator 71 is given by 

(3.22) A = -V 3 a^^a 2 ^)^ - V 2 a 2 ^)^ _ u 2 a^^a^)^ - Ife ofx)^, 
where we have defined four group parameters Q 



V3 — 



2 

Inserting (|3.21l) into (|3.17l) we find 

(3.23) Au afl = {-d t + {L 2 )) Ul . . 

Given an expression for uo,o and appropriate BCs, one can use PDE ()3.23j) to find an expression for ui,q. 
This is as far as we will take the analysis of equation (|3.8j) . We now return to the O(VS) equation (|3.9j) . 

Singular Perturbation Analysis of Equation (|3.9[) 

The singular perturbation analysis of (|3.9[) proceeds by inserting expansions (|3.6[) and (|3.7[) into (|3.9p and 
collecting terms of like-powers of y/e. The resulting order 0(y/S/e) and 0(VS/y/e) equations are 

(3.24) 0{VS/e): = £ «o,i, 

(3.25) Q(y/8/yfe) : 

where we have used M3Uq,o = 0. We note that if uo,i and are independent of y, equations (|3 . 24[) and 
()3.25j) will automatically be satisfied. Thus, we choose uo,i = uo,i(x, z) and u\^\ = ui t i(x,z). Continuing 
the asymptotic analysis, the order 0(VS) equation is 

(3.26) 0(VS): = £ ii 2 ,i + (-9 t +£ 2 )wo,i+MiMo,o, 



4 The phrase group parameter refers to any z-dependent parameter which can be calculated as a moment of model-specific 
functions. As we shall see, the effect that the functions (/3, A, g, T) have on the approximate price of a derivative asset is felt 
only through eight group parameters. 
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where we have used £1^1,1 = and Matt^o — 0. We note that equation (|3.26p is a Poisson equation for W2.1 
of form (|3.14p . By the Fredholm alternative, in order for p.26[) to admit a solution ui,\ € L 2 (M.,tt) centering 
condition (|3.15[) must be satisfied. In p.26[) centering condition (13. 15)) corresponds to 

(3.27) = {-d t + (L 2 )) uo,i + (Mi) u 0fi . 

Note that Uo,o(*) x i z ) depends on z only through a(z) and fil(z). Thus, in (|3.27l) (Mi) can be written 

(Mi) = -m z , 

(3.28) B = -Viaf*)^ - V 0l Vi := g Pxz (/) , V := -9 (V) , 

(3.29) d, = o , d ir + ~0dj t[ , a';=d z a, fit := d z ~JU. 

Note that we have introduced four more group parameters: Vi, Vq, a 1 and . This is as far as we will 
take the asymptotic analysis of equation (|3.1[) . For convenience, we review the most important results of 
this section. 



Main Results of the Asymptotic Analysis 

(3.30) 0(1): (-d t + (£ 2 )) Uo ,o = 0, u ,o(0, x, z) = H(x), 

(3.31) 0(Vi): (-ft + <-C2))ui,o=^uo,o, Ui j0 (0,a:,*) =0 

(3.32) Q(VS): {-d t + (Z 2 )) «o,i = 30»u Ol o, u ,i(0,x,z) = 0. 

The operators (£ 2 ), A 2 and d 2 are defined in (|3TT5|) . ([3~2"2"|) . p~2g|) and p~29"j) respectively. Note that we 
have imposed BCs at t = 0. 



3.2 Explicit Solutions for uq.o, u\fl and «o,i 

In this section we shall explicitly solve equations (|3.30[) . f|3.31[) and (|3.32l) in terms of the eigenfunctions {ipn} 
and eigenvalues {A„} of the operator (£2)- To begin, we note that (£2), given by (|3. 181) . has the form of an 
infinitesimal generator of a one-dimensional diffusion (jl.3l) with volatility aa(x), drift (b(x) — fil a(x)) and 
killing rate k(x). The dom((£2)) includes BCs, which must be imposed at the endpoints ei and e 2 . Appendix 
IA.2I describes the appropriate BCs to impose for a general one-dimensional diffusion with a generator of the 
form (TQ1) . 

Throughout this section we assume (£2) has a purely discrete spectrum. We fix a Hilbert space "K = 
L 2 (I,m) where m is the speed density corresponding to (£ 2 ). The operator (£ 2 ) is self-adjoint in Jf and its 
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domain is a dense subset of !K. Thus, the eigenfunctions {ip n } of (£2) form an orthonormal basis in JC. It 
is not necessarily true that either A : JC — > JC, 23 : 5£ — > "K or 9^ : "K — > "K. As such, we define 

dom(yi) := {ip £ 3£ : Atp £ "K\ , dom(S) :={)/)£?{: 23^ G , 

dom(<9 z ) := {V> S % : d z ip E %] . 

Theorem 3.1. Assume that we can solve the following eigenvalue eguation 

(3.33) - (£ 2 ) V>n = K.ipn, 4>n e dom((£ 2 )) , 

and assume H e M. TTien </ie solution i*o,o to (|3.30[) is given by 

UQfl = ^ Cn1p n T n , Cn = (lp n ,H), T„ = e~* A ™. 

Proof. One can easily verify that uo,o satisfies PDE (|3 . 30[) assuming (|3.33l) holds. To see that the BC 
ito.o(0, x, z) = H(x) is satisfied, notice that T n (0, z) — 1 and apply equation (| A.2|) to the payoff function H 

Id H = 2J (lpn,H) 1p n = ^ °n^n- 
n n 

□ 



Theorem 3.2. Let c„, if> n a71 ^ ^ e as described in Theorem \3.1\ and define 

A k , n ■= (ipk,Ai/) n ) , U k ,n ■= t~ — 7^- 

T/ien i/ie solution «i,o to equation (|3.31[) is 
(3.34) iti,o = X X 

n k^n n 

Proof. See appendix IA. 51 □ 
Note that ui t o is linear in the group parameters (V3, V2, U2, Hi). 

Theorem 3.3. Let c n , ip n and T n be as described in Theorem \3.1l let Uk, n be as described in Theorem 
and define 

S fcl „:=(Vfc,20^n), 2 M :=(fc%), ^ fc|W := , T *~ T " 2 + , * T " ■ 

(Aft - A„) A fc - A„ 

Then the solution 1*0,1 to equation (|3.31[) is 



"0,1 = C n' S 'k,n'>PkUk,n ~ X C ™ <T ™ 



n fc^n 



(3.35) + X X Cri^fe.n'/'fe (<9 Z A„) Ffc,„ + Cn'Bn^n (d z X n ) \t 2 T n . 

n k^n n 
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Proof. See appendix I A. 61 □ 
Note that uo,i is linear in (Via', VifSl , V ct', VqJU'). 

Accuracy of the Pricing Approximation 

We have now derived an approximation u e ' S « 2/0,0 + \/^ u i.o + \^ u o,i for the price of a derivative- asset. 
However, this derivation relied on formal singular and regular perturbation arguments. In what follows, we 
establish the accuracy of our approximation. For our accuracy result, in addition to the assumptions listed 
in section [A. 3| we shall need one additional assumption 

• The payoff function H(x) and all of its derivatives are smooth and bounded. 

Obviously, many common derivatives - e.g., call and put options - do not fit this assumption. To prove the 
accurac y of our pricing approximation for calls and pu ts would require regularizing the option payoff as is 



done in 



Fouque. Papanicolaou. Sircar, and Spina 



(|2003l ). The regularization procedure is beyond the scope 



of this paper. As such, we limit our analysis to options with smooth and bounded payoffs. Our accuracy 
result is as follows: 

Theorem 3.4. For fixed (t,x,y,z), there exists a constant C such that for any e < 1, S < 1 we have 

u e ' 5 - f«o,o + y/eu Qt0 + Vduo t ij < C (e + S) . 
Proof. See appendix IA. 71 □ 

Theorem l3.4l gives us information about how our pricing approximation behaves as e — > and 5 — > 0. In 
practice, both e and 5 are small, but fixed (they do not go to zero). Without knowing what the constant 
C is in theorem 13.4) it is difficult to gauge exactly how good our pricing approximation is. As such, in the 
examples provided in section |4j we will compare the approximate prices of derivative- assets (calculated using 
Theorems 13.11 13.21 and I3.3[) to their exact prices (calculated via Monte Carlo simulation) . 



4 Examples 

In this section we compute the approximate price of three derivative-assets: a double-barrier call option, a 
bond in a short-rate model, and a European call on a defaultable stock. 
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4.1 Double-Barrier Call Option with Multiscale Stochastic Volatility 

In our first example, we let X represent the value of a non-dividend paying asset (e.g., a stock, index, etc.). 
Often, X is modeled as a GBM with constant volatility (e.g., Black-Scholes). Here, we model X as a GBM 
with multiscale stochastic volatility. Specifically, the P dynamics of X are given by 



dX t = rX t dt + f (Y t , Z t ) X t dW t x , 



h(X t ) = 0, 



where r is the risk-free rate of interest and Y and Z are fast- and slow-varying factors of volatility, as 
described in (12. 4|) . Note that, as it should be, the discounted price of the asset (e~ rt X t ) is a martingale 
under P. We will calculate the approximate price of a double-barrier call option written on X. 

To start, we use equations (|1 .4[) and (|3 . 1 8[) to write the operator (£2) and its associated speed density 
m(x) 



(4.1) 



(£2) = \f 2 x 2 d 2 xx + rxd x -r, 



m { x ) = =2 — 2 eX P 1 =2 l0 S X 



For a double-barrier call option with knock-out barriers at L and R, the option payoff is 



H(X t )I {T>t} = (X t -K) + l 



{Tj>t}> 



I=(L,R), 



< L < K < R, 



To calculate the value of this option we must first solve eigenvalue equation (|3.33[) with (£2) given by (|4.1[) 
and with BCs 



lim tpn(x) = 0, 



lim ip n (x) = 0. 

x/*R 



Note that we have imposed the regular kill ing BC at the e ndpoints L and R. The solution to (|3.33p with 



the above BCs can be found on page 262 of 



Linetskvl (j2007l ) 

n^/x f-r \ . fmr\og(x/L) 

ip n (x) = — z exp — o logx sin — — ; — 7—— 



A n — — 



n = 1,2,3, 



r a 

v ~ i ~ 2' 



2 Viog(i?/i) / ' V 2 

Next, we use expressions (|3.22[) and (I3.28|) to write expressions for the operators A and 23 



(4.2) 



A = -V 3 xd x x 2 d 2 xx -V 2 x 2 d 2 xx , 



V = -V 1 xd x -V Q . 
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Using (|4.2[) it is now straightforward to calculate inner products Ak, n , 23fe,n and "Bk,n- For k 7^ n we find 

(-1 + (-l) k+n ) fcn (4n 2 7T 2 ff 4 + (-12r 2 + 4m 2 + ct 4 ) log 2 (i?/L)) \ 



23fc,n = Vl 



2(fc 2 -n 2 ) cr 4 log 3 (R/L) 

4(-l + (-l) fe+n ) knr\ 
(fc 2 -n 2 )<j 2 log(i?/L) J ' 

2 (-1 + (~l) k+n ) kn 



V 



•B k ,n = -Via' (T 



(fc-n)(fe + n)log(ii/i)' 

, /8 (-1 + (-l) fc+n ) knrlog(R/L)\ 



Vo<? 



T fc , 



(fc 2 -n 2 ) 2 7r 2 CT 3 
4fcnr (log(L) - (-l) fc +™ log(i?)) 



(fc 2 -n 2 )a 3 \og(R/L) 
2 (-1 + (-l) fe+n ) fcn ((fc - n)(k + n^W 4 - 2r (-2r + a 2 ) log 2 (i?/L)) 
(fc 2 - n 2 ) 2 tt 2 ct 5 \og{R/L) 



and for k — n we find 



A-n.n — — V3 



1 / 3n 2 7r 2 ^ 



3 =3 



a 3 Vlog'WL) 



— 9 

(7 



9 



log 2 (i?/L), 



- V, 



a 2 



log* (R/L) J a 



K 



23«,n = Vi 



2r- a 2 
2a 2 



Vo, 



S ™ (\o£(R)-W{L)) 

n ' n ~ 1(7 \2<f W 4 log (R/L) j 

The calculation of c n can be found on page 262 of 
c n = (M-),(--K) + ) 



_, f 1 r(log 2 (E)-log 2 (L)) 
U + a 3 log(i?/i) 



Linetskvl (|2007D 



w 2 + 



log (i2/L) 

exp (&z) (cu n cos (w„.ft) — z sin (w„.S) ) — exp (ilz) (— l)"o; n J , 



1 



w« := -TT-, i?:=~log( — ), U:=-log( v 



K 



R 



Approximate option prices can now be computed using Theorems 13. 1[ 13.21 and 13.31 

On the left side of figure Q] we plot the approximate price uo,o + v^ u i.o °f a double-barrier call option 
for a specific model that has only a fast-varying factor of volatility. We suppose the dynamics of Y and the 
volatility function / are given by 

cr exp (It) 



(4.3) 



dY t = ( --Y t - -LjSErf(y t ) ] dt + fidWi 



y 

t > 



cxp(-/? 2 /2)' 



Er%) 



e"' dt. 



From comparison we also plot the full price u e (calculated by Monte Carlo simulation) and Uo,o> which 
corresponds to the Black-Scholes price with volatility cf. On the right side of figure [T] we plot the approximate 
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price uo,o + V^uo,i of a double-barrier call option for a specific model that contains only a slow- varying factor 
of volatility. We suppose the dynamics of Z and the volatility function / are given by 



(4.4) 



dZ t = (-8Z t -V6g Eri(Z t )) dt + g dW t 



f(Zt) 



a exp(Z t ) 



exp (z) 

For comparison, we also plot the full price u 5 (calculated by Monte Carlo simulation) and the Black-Scholes 
price uo,o- As expected, as e and 6 go to zero, the approximate price converges to the full price, which 
conveges to the Black-Scholes price. 



4.2 Vasicek Short-Rate with Multiscale Stochastic Volatility 

In our second exampl e, we let X rep resent the short-rate of interest. One of the most widely known short- 



rate models is that of 



Vasicekl (1977) , in which X is modeled as an OU process. Here, we model X as an 



OU with multiscale stochastic volatility. Specifically, the P dynamics of X are given by 

dX t = ^ K (e-X t )-f(Y t ,Z t )n(Y t ,Z t )yt + f(Y t ,Z t )dW t x , r(X t )=X t , h(X t ) 



0, 



where Y and Z are fast- and slow- varying factors of volatility, as described in (|2.4p . We will calculate the 
approximate price of zero-coupon bond in this setting. [f| 

To start, we use equations (| 1 .4[) and (|3.18[) to write the operator (£ 2 ) and its associated speed density 
m(x) 



1 



2-expl 
a \ a 



(4.5) (L 2 ) = -a l di x +K{9-x)d x -x, m(x) 
For a zero-coupon bond, the payoff at maturity is 

(4.6) H(X t )l {T>t} = l. 

In order to price a bond with payoff (|4.6[) , we must solve eigenvalue equation p.33[) on the interval 
I = (—00,00) with (£2) given by (|4.5[) . As both —00 and 00 are natural boundari es, no BCs need to be 



specifi ed. The solution to this eigenvalue problem can be found in equation (4.6) of 
I2OO4J) 



Gorovoi and Linetskv 



i> n = K n exp [-AS,- -A 2 ) H n (f + A) , 



7T 2" +1 n! 



1/2 



f 3/2 ' 



An — An 



2k? 



a x 



n = 0,1,2, 



We note that r(Xt) = Xt may become negative whe n X is described by an OU pr ocess. As such, one may wish to impose a 



reflecting boundary condition at x = 0, as carried out in 



Gorovoi and Linetskv 



12004). However, as an OU without a reflecting 



boundary is most prevalent in literature, this is the case we treat here. 
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Here, {H n } are the (physicists') Hermite polynomials. Next, we use (|3.22p and p. 281) to write expressions 
for the operators A and CB 



A = -Va dl xx - (V 2 + U 2 ) d 2 xx - Ui d x , 



'B = -V 1 d x -V . 



It is now straightforward to calculate inner products Ak, n , ^k,n and T$k,n- Using the recursion relations 



d x H n = 2 nH n -i, 



2 x H n = H n+ i + d x H n , 



we find 



A k 



3Am 



MY, 



m=0 \ m 



■n 3 - m ciy/a 



(n — to)! N„_ 



-<5fe.n-r 



2An / 2 

(V 2 + U 2 )<|^ 



m=0 \ TO 



-l\ 2 ~ m /2^ 



Hi 



T$k,n — — Vl ( — ) £fc,n 



-Vict' 



«/ V * / (n-l)!^-! 

a ) (n-l)!^/^" 1 

-l\ / 1 71 

k J \2a k 3 a 

zl\ [AS] + ( hlE\ (L-^-l 

K J \ n 3 / 2 / V a \2a k 3 a 



(n — m)! 'Nn- 
5k,n-l\ > 

Vq <5/c,n, 



-5k,n-r 



5k,n 



+ 
+ 

V a' 
+ 

v{Jn' 



\ G 

a J { k- - 



K ) ( a 

\ a )\ a )\ (n-3)!K„_ 3 
1 a n\ , ( 4 
2cf k 3 
2\ n\K 



(n - 2)!N„ 



(n-l)!X„_! 

<5fe,n-2 



= I <Sfc,n 



K 3/ 2 y ( n _ i)!N„_ x 



a / (n-2)!K„_ 2 



4 



fe,n-2 



1 

- I 

K 2 



4,™ + 



5k,' 



a^ 3 / 2 ; (n-l)!X n _i 
2 \ n\N„ 



a 2 J (n-2)W n -2 



-5k,n-2 



o^k) (n - l)!3\f n _i 
The computation of c n be found on page 63 of in 



5 k,; 



Gorovoi and Linetskvl (|2004l ) 



Cn=W„,l) = =«AV^ 4 . 

a v k 

The approximate price of a bond can now be calculated using Theorems 13.1! 13.21 and 13.3 
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Yield Curve 

For a zero-coupon bond, it is often the yield curve, rather than the bond price itself, that is of fundamental 
importance. The yield R e ' S of a zero-coupon bond that pays one dollar at time t is defined via the relation 

u c ' s = exp (-R e > s t) . 

We can obtain an approximation for the yield of a zero-coupon bond by expanding both the bond price u e ' S 
and yield R e ' S in powers of y/e and y/6 as follows 

u ,o + V~em, + Vduo,! + ■■■ = e -Ko+v^ 1 ,oW5iJo, I +-) t 

= e - R "' ot + v/e (-Ri,o t) e - Ro ^ + VS (-Ro,i t) e - Ro - at + ■■■ . 

Matching terms of like-powers of i/e and V~5 we obtain 

i? £ ' 5 « Ro,o + VeRi.o + VSRo,u 

D 1 1 I \ D ~ M 1,0 D — "0,1 

^io. o — — 7 log (U q) , A1,0 = ) tto.i = , 

tu 0fi tu 0fi 

On the left side of figure [2] we plot the approximate yield i?o,o + V^-Ri.o of a zero coupon bond for 
a specific model that has only a fast-varying factor of volatility. We suppose the dynamics of Y and the 
volatility function / are given by (|4.3j) . For comparison, we also plot the full yield R e (calculated by Monte 
Carlo simulation) and the Vasicek yield Rq.o- O n the right side of figure [5] we plot the approximate yield 
^o,o + Vfi Ro.i of a zero coupon bond for a specific model that has only a slow- varying factor of volatility. 
We suppose the dynamics of Z and the volatility function / are given by (|4.4[) . For comparison, we also plot 
the full yield R s (calculated by Monte Carlo simulation) and the Vasicek yield i?o,o- As expected, as e and 
S go to zero, the approximate yield converges to the full yield, which converges to the Vasicek yield. 

4.3 Jump to Default CEV with Multiscale Stochastic Volatility 

In our final example, we consider a non-dividend-paying, defaultable asset St — I{ T>t }X t . As S must be 
non-negative, we let the state space of X be (ei,e2) = (0, oo). We base our multisc ale diffusion on the 
jump to default constant elastic variance model (JDCEV) of 
dynamics of X prior to default are given by 



Carr and Linetskvl (|2006l) . Specifically, the 



dX t = 



(/x + c V t 2 ") X t dt + (f(Y t , Z t ) A7) X t dW f x , h(X t ) =f i + cX 2 t r >. 



For computational convenience we have set the risk-free interest rate to zero: r = 0. The constants \x and c 
are assumed to be strictly positive. As always, Y and Z are fast- and slow- varying factors of volatility, as 
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described in (|2.4[) . Note that the volatility of X has both a local component and a nonlocal multiscale 
component f(Y tl Z t ). We assume n < so that the local component of volatility X% increases as X t decreases, 
reflecting the fact that price and volatility are negatively correlated. The stochastic hazard rate h(X t ) also 
increases as X decreases, capturing the idea that the probability of default increases as X tends to zero. 
Note that S is a P-martingale, as it should be. We will calculate the approximate price of a European put 
option written on S. The price of a European call option can be obtained through put-call parity. 
To begin, we use (|1.4[) and f|3.18[) to write the operator (£ 2 ) and its associated speed density m(x) 

(4.7) (L 2 ) = ^a 2 x 2r ' +2 d 2 xx + + cx 2 ^)xd x - (/x + cx 2 ^) , 

m ( x ) = i s x 2c / w2 - 2 - 2r iexv{Ax- 2r >), A= ^ 



J 2 n " <t 2 M 
For the diffusion associated with infinitesimal generator (|4.7p the endpoint e 2 = oo is a natural boundary. 

However, the classification of endpoint e\ — depends on the values of n and c/a 2 . The classification is as 
follows 

c/a 2 > 1/2 and i] < 0, e\ = is natural, 

cja 2 e (0, 1/2) and n € [c/ct 2 - 1/2, 0), ei = is exit, 

cja 2 € (0, 1/2) and 77 < c/ct 2 - 1/2, ei = is regular. 

If the parameters (c, a, 77) are chosen such that e\ = is regular, then we specify e\ = as a killing boundary. 
To calculate the approximate price of a European put we must solve the eigenvalue equation (|3.33j) on the 
interval (0, 00) with (£2) given by (|4.7|) and with the BC 

lim ifj n (x) = 0, if c/a 2 G (0, 1/2). 

x\0 



The solution is given in equation (8.11) of Theorem 8.2 in iMendoza-Arriaga. Carr. and Linetskvl (|2010l ) 



i>n = ^ /2 y ^T^ X 6XP (~ A ^ 2 ") 4^1 (^^ 2 '0 > n = 1,2,3,-. , 

1 + 2 {c/a 2 ) 



where {Ln } are the generalized Laguerre polynomials. Next, we use (|3.22|) and p.28p to write expressions 
for the operators A and 23 

^ = -V 3 x n+1 d x x 2v+2 d 2 - V 2 x 2ri+2 d 2 V = -V! x r,+1 d x - V . 



21 



Analytic expressions for Ak, n , ^>k,n and H>k,n are easily derived by making the change of variables Ax 2r) — > y, 
using d y L v n {y) = -L [ ^\y) and 

y^e-yL^\y)L ( i\y)dy 

_ r(a-7 + n)r( / 8 + l + m)r(7 + l) / 7 + 1, -m, 7 + l-a 



r(a-7)r(^ + l)n!m! 



(3 + 1, 7 + 1 — a — n, 



where whe r e v F a is a generalized hypergeometric function (the above formula is given in equation (14) of 



Shawagfehl (|201lh ). As the formulas for Ak, n , 25fe,n and T>k,n are quite long, for the sake of brevity, we do 
not provide them here. 

The payoff of a European put option with strike price K > can be decomposed as follows 
(4.8) (K - S t )+ = {K- X t )+ l {T>t} +K(l- l {T>t} ) . 

The first term on the RHS of (|4.8[) represents the payoff of a put given no default prior to time t. The second 
term represents the payoff of a put option given a default occurs prior to time t. Thus, the value of a put 
option with strike price K — denoted u e ' S (t, x; K) - can be expressed as the sum of two parts 

u e ' 5 (t, x; K) = u Q ' S (t, x; K) + u% 5 (t, x; K), 

where 

u e ' S (t, x- K) - E x>y>z [{K - X t )+l {T>t} ] , 



e,<5 



{t,x;K) = K-KE x>y , z [l {r>t} ] 

= K-K % XiViZ [6 X , (X t ) l {T>t} ] dx' 
Jo 



(4.9) =K-K I u{\t,x-x')dx\ 







e,<5 

U 



1 (t,x]x') =E X ,„,2 [5 X > (X t ) I{ T >t}] • 

Note, because 1 ^ L 2 (M + ,m), we have used the fact that 1 = f^ 8 x i{X t ) dx' on the set {r > t}. This 
substitution comes at a cost; the integral in (|4.9[) must be computed numerically. However, numerical 
evaluation of (|4.9[) is not computationally intensive and does not pose any major difficulties. 

Since the payoff functions H (x) = (K — x) + and Hi(x) = 5 x >(x) belong to L 2 (R + , m), we may calculate 

Co,n = (lpn{-), (k - + ) i c l,n = (^n, 4') • 
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The ex pression for co. n can be found in equation (8.15) of Theorem 8.4 in lMendoza-Arriaga. Carr. and Linetskv 



(|20 10T ) . The expression for ci j7l is computed trivially. We have 
A»l 2 + l K 2c l° 2+1 - 2 ^T{ V + n) 



c 0,n — 




r(v + i)(n-i)i L(: + 1)(AK _ 2v) 



{c/a*) + \r,r* \ „ + L £^! +2 / Tw + n :i) 

Cl,n = 4>n(x')m{x'). 

The approximate price of a European put option can now be computed using Theorems 13. 1[ 13.21 and [ 

For European options, it is often the implied volatility induced by an option price, rather than the option 
price itself that is of primary interest. Recall that the implied volatility I e ' S of a put option with price 
u e ' S (t, x; K) is defined implicitly through 

u e ' S (t,x;K) = u BS {t,x,I e ' S ;K) 

where u BS (t,x, I e S ; K) is the Black-Scholes price of a put as calculated with volatility J e ' . 

On the left side of figure [3]we plot the implied volatility induced by the approximate price «q,o + y/euifi 
of a put option for a specific model that has only a fast-varying factor of volatility. We suppose the dynamics 
of Y and the volatility function / are given by (|4.3p . For comparison, we also plot the implied volatility 
induced by the full price u e (calculated by Monte Carlo simulation) and the implied volatility induced by the 
JDCEV price ito.o- O n the right side of figure [3] we plot the implied volatility induced by the approximate 
price ito.o + V^uo.i of a put option for a specific model that has only a slow- varying factor of volatility. We 
suppose the dynamics of Z and the volatility function / are given by (|4.4[) . For comparison, we also plot 
the implied volatility induced by the full price u s (calculated by Monte Carlo simulation) and the implied 
volatility induced by the JDCEV price 1*0,0- As expected, as e and S go to zero, the implied volatility induced 
by the approximate price converges to the implied volatility induced by the full price, which converges to 
the implied volatility induced by the JDCEV price. 



5 Review and Conclusions 

This paper develops a general method for obtaining the approximate price for a large class of derivative- 
assets. The payoff of the derivatives may be path-dependent and the process underlying the derivative-assets 
may exhibit jump to default as well as combined local/nonlocal stochastic volatility. The intensity of the 
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jump to default event may be state-dependent and the nonlocal component of volatility may be multiscale, 
driven by one fast- varying and one slow- varying factor. 

One key advantage of our pricing methodology is that, by combining techniques from spectral theory, 
singular perturbation theory and regular perturbation theory, we reduce the derivative pricing problem 
to that of solving a single eigenvalue equation. Once this equation is solved, the approximate price of a 
derivative-asset may be calculated formulaically. We have illustrated the simplicity and flexibility of our 
method by calculating the approximate prices of thre derivative assets: a double-barrier option on a non- 
defaultable stock, a European option on a defaultable stock, and a non-defaultable bond in a short-rate 
model. 

We believe that the flexibility of our framework, as well as the analytic tractability that our pricing 
methodology provides merit further research in this area. A logical next step, for example, would be to 
extend the results of this paper to include cases where the eigenvalue equation (|3.33p does not have a purely 
discrete spectrum. 
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A Appendix 



A.l Self- Adjoint Operators acting on a Hilbert Space 



In this appendix we summarize some basic properties of self-adjoi nt operators acting on a Hilbert space. A 



detail ed exposit i on on this topic (including proofs) can be found in 



follow 



Reed and Simonl (1980). We shall closely 



Linetskvl (|2007j ). who provides a more streamlined review. 



Let !K be a real, separable [j Hilbert space with inner product (•, •). A linear operator is a pair (dom(£), L) 
where dom(£) is a linear subset of !K and L is a linear map L : dom(£) — > J£. The adjoint of an operator 
L is an operator £* such that (£>f,g) = (/,£*#), V/ <E dom(£),p G dom(£*), where 



dom(£*) := {g e JC : 3h e 5£ such that (Zf,g) = (f,h) V/ G dom(£)}. 



3 A Hilbert space is separable if and only if it admits a countable orthonormal basis (i.e., Schauder basis). 
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An operator (dom(£),£) is said to be self-adjoint in J{ if 



dom(£) = dom(£*), 



(Lf,g) = (f,Lg) 



V/,5 G dom(£). 



Throughout this appendix, for any self-adjoint operator L, we will assume that dom(£) is a dense subset of 



Given a linear operator £,, the resolvent set p(£>) is defined as the set of A G C such that the mapping 
(£ — Id A) is one-to-one and R\ := (£ — Id A) -1 is continuous with dom(i?A) = The operator R\ : "K — > 5£ 
is called the resolvent. The spectrum o~(£) of an operator £ is defined as o~(L) := C\/?(£). If L is self-adjoint, 
its spectrum is non-empty and real. We say that A G o-(L) is an eigenvalue of L if there exists "0 S dom(£) 
such that the eigenvalue equation is satisfied 



A function ^ that solves (|A.1[) is called an eigenfunction of L corresponding to A. The multiplicity of 
an eigenvalue A is the number of linearly independent eigenfunctions for which equation (jA.lj) is satisfied. 
The spectrum of an operator L can be decomposed into two disjoint sets called the discrete and essential 
spectrum <r(£) = o~d{&) U o~ e (L). For a self-adjoint operator £, a number A G M belongs to crd(£) if and only 
if A is an isolated point of er(£) and A is an eigenvalue of finite multiplicity. 

The spectral representation Theorem is an important tool for analysing self-adjoint operators acting on 
a Hilbert space. We state this theorem below in a form which is convenient for the computations in this 
paper. 

Theorem A.l. Assume L is a self-adjoint operator in IK and assume L has a purely discrete spectrum (i.e., 
o~ e (L) = {$})■ The Spectral Representation Theorem states that Lf has an eigenfunction expansion 



where the sum runs over all solutions {\ n ,ipn} of the eigenvalue equation (|A.1[) . Furthermore, for any 
real-valued B or el- measurable function on M one can define an operator </>(£) using functional calculus 



(A.l) 



L i\) = Xijj. 




V / G dom(L) 



n 




V/Gdom(0(£)), 



n 



dom(0(£)) :={/Gdom(£) :^0 2 (A n )(^„,/) 2 <oo}. 



n 



The operator is self-adjoint in !H and dom(0(£)) C dom(£). 



Proof. See lReed and Simonl (|l980l ) Theorem VIII.6. 



□ 
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Note that setting 0(A) = Id yields 

(A.2) Id/ = X>n,/)V>n, 



which is equivalent to saying that the eigenfunctions {tp n } of a densely defined self-adjoint operator in "K 
form a Schauder basis. In fact, the basis can be chosen to be orthonormal (ipmipm) = S n . rn . Also note, 
setting (j)(X) = R\ yields an eigenfunction representation of the resolvent operator 



A„-A 

A.2 Boundary Conditions 



4>n 



V/eM,Ae p(L). 



According to 



Fellerl (11954) . the endpoints e\ and e 2 of a one-dimensional diffusion in an interval I can be classi- 



fied as either natural, ex it, entrance or regular. The classification, which can be found in 



Borodin and Salminen 



(|2002l ); iLinetskvl ([20071 ). is done as follows. For a general infinitesimal generator L of the form (II .3[) one can 

associate a scale density 



s(i): == p (-£§l4 



(scale density) 



where the lower limit of integration xq € (ei,e 2 ) may be chosen arbitrarily. From s one can define a scale 
function § 

rV 



S ([*,»]) 



s(z) dz, 



x,y e (ei,e 2 ) 



S([x,e 2 )) := lim §([x,y]) 



Note that the above limits may be infinite. For some arbitrary y £ (ex, e 2 ) we define 



h 
Ji 



S ((ei,x]) (1 + A:(a;)) m(a;) dx, 



S ([a;, y]) (1 + k(x)) m(x) dx, 



h 
■h 



§ ([x, e 2 )) (1 + k(xj) m(x) dx, 



S ([y, x]) (1 + k(x)) m(x) dx. 



An endpoint is classified as 

• Natural if 7j = oo and Ji — oo. No BC needs to be specified at a natural boundary. The interval I is 
taken to be open at a natural boundary. 

• Exit if Ii < oo and Ji = oo. The appropriate BC at an exit boundary is 

lim ij){x) = 0. 

x— yei 

The interval / is taken to be open at an exit boundary. 
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• Entrance if Zj = oo and Ji < oo. The appropriate BC at an entrance boundary is 

lim = o. 

The interval 7 is taken to be open at an entrance boundary. 

• Regular if Zj < oo and Ji < oo. We must specify the behavior of a diffusion at a regular boundary. 
Here, we consider only killing and instantaneously reflecting behavior, for which the appropriate BCs 
are 

lim ip( x ) = (killing BC), lim x — (instantaneously reflecting BC) 

x— let x^fCi S[x) 

The interval I is taken to be open at a regular boundary specified as a killing boundary and closed at 
a regular boundary specified as instantaneously reflecting. 

The domain of L is then 

(A.3) Dom (L) = {/ e L 2 (I, m) : /, 8 x f e AC loc (I), Lf e Z 2 (I, m) , BCs at d and e 2 } , 

whe re ACi oc (I) is t he space of functions that are absolutely continuous over each compact subinterval of I 



see 



Linetskvl (|2007h . p. 242). The BCs at e\ and e 2 correspond to the BCs specified above for natural, exit, 



entrance and regular boundaries. 

A.3 Specific Model Assumptions 

1. We assume existence and uniqueness of (X, Y, Z) as the strong solution to (|2.ip . 

2. We assume existence and uniqueness of (X, Y, Z) as the strong solution to (|2.4j) . 

3. There exist positive constants Ca < oo and Cr < oo such that HAHoo < Ca and ||r||oo < Cr- 

4. Define the time-rescaled process Y t (1) :=Y et . Under P, the process Y^ has infinitesimal generator fo- 
under P we assume: 

(a) The process Y^ is ergodic and has a unique invariant distribution n with density n. 

(b) The operator Lq has a strictly positive spectral gap - meaning the smallest non-zero eigenvalue 
A m in of (—Lq) is strictly positive. 

(c) The process Y^ is reversible - meaning Lq is self-adjoint acting on L 2 (M, tt). 
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These assumptions guarantee (see lFouaue. Papanicolaou. Sircar, and Solnal (|201ll ). p. 93) exponential 
convergence of Y^ 1 ' to its invariant distribution 



E 



9 Y t 



(i) 



(g) < C exp {-\ min t) . 



Vg £ L 2 (K,tt). 



The above assumptions also ensure (see lFouaue. Papanicolaou. Sircar, and Solnal (|201ll ), p. 139) that 



for all k £ N there exists C(k) < oo such that 



(A.4) 



supE 



Y, 



(i) 



< C(k). 



5. Define the time-rescaled process Z t := Z t /g. Under P, the process Z^ has infinitesimal generator 
M2. Under P we assume the process Z^ 1 ' admits moments that are uniformly bounded in s < t. That 
is, for all k £ N there exists C(t, k) < 00 such that 



(A.5) 



supE 

s<t 



zP 



< C(t,k). 



6. We assume that the functions f(y,z) and Q,{y,z) satisfy <J 2 (z) < 00, fil(z) < 00 and the solutions 
4>(y, z) and rj(y, z) to Poisson equations (|3.19[) are at most polynomially growing. 

7. The functions a(x) b(x), r(x) and h(x) satisfy a > 0, a £ C 2 (I), b £ C 1 (/), r > 0, r £ C(I), h > 0, 
and h £ C{I). 

8. The spectrum of the operator (£2), defined in (13. 18[) . is simple and purely discrete. 

We note that two of the processes that are most commonly used to model volatility - the Cox-Ingersoll-Ross 
(CIR) and Ornstein-Uhlenbeck (OU) processes - satisfy the assumptions placed on both Y^ and Z^ . 



A.4 Poisson Equations and the Fredholm Alternative 

In this appendix we review the existence and uniqueness of solutions to Poisson equations. Central to 
this discussion will be a statement of the F redholm alternative. Our presentation follows page 93 of 



20071 ). 



Fouaue. Papanicolaou. Sircar, and Solnal (|201ll ). as well as page 124 of lFouaue. Gamier. Papanicolaou, and Spina 



Let £ be a self-adjoint operator densely defined on a real separable Hilbert space Di, and let {4>n, A„} be 
the complete set of solutions to eigenvalue equation £?p n = X n ip n . Consider the following Poisson problem: 
find, i\) £ r K such that 

(A.6) (£-A)V = X, 

where the function x € 3~C and the constant A are given. 
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Theorem A. 2. The Fredholm Alternative states that one of the following is true: 

1. Either A is not an eigenvalue of L, in which case equation (|A.6I) has a unique solution 



2. Or, A is an eigenvalue of L. Suppose this is the case. Let A = Ai = A2 = • • • = A m (i.e., the eigenvalue 
A has multiplicity m). Then (|A.6|) has a solution if and only if (jp n ,x) = f or a ^ n < m. Assuming 
(VVu x) = f or a M n < m, a solution to (|A.6I) has the form 

n>m n<m 



Proof. See 



Reed and Simonl (|1980h . Theorem VI. 14 and the ensuing corollary. □ 



Classically, the Fredholm alternative Theorem holds for compact operators on a Hilbert space. However 
the Theorem also holds true for differential operators £ of the form (|1.3|) , with domain (|A.3j) ac ting on the 



Hilber t space "K = L 2 (I,m), where m is the speed density corresponding to L (see section 9.4.2 of 



Haberman 



l|2004h ). 

In particular, we note that A = is an eigenvalue of L , which is a self-adjoint operator in L 2 (R,7r). 
The corresponding (normalized) eigenfunction is the constant ip\ — 1. Thus, in order for Cqu = % to have a 
solution u £ L 2 (R,7r) we must have (l,x) = / X^dy =: (x) = 0, which is the centering condition (|3.15|) . 

A . 5 Proof of Theorem ET21 

We must show that Ui,0) given by (|3.34p satisfies PDE and BC (|3.3ip . It is obvious that the BC ui.o(0, x, z) = 
is satisfied. To show that u%fi satisfies PDE (|3.3ip we note that 



(A.7) Au 0fi = ^2 c n (Aip n ) T n = ^2 ^2 c n A k:Tl ip k T n , 

n n k 

where we have used (|A.2[) in the second equality. Now, using (|3.33p and 

(A.8) (-d t - A fc ) XJ Kn = T„, (-d t - \ n )tT n = -T n 

it is easy to show that 

(-d t + (L 2 ))ui, = CO}. 
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A.6 Proof of Theorem [3731 



We must show that uq,i, given by (|3.35l) satisfies PDE and BC (|3.32l) . It is obvious that the BC u Oj i(0, x, z) 
is satisfied. To show that uo,i satisfies PDE (|3.31[) we note that 

C n (m z ^ n ) T n + ^ (d z C n ) (S^ B ) T n + ^ C„ (S^ n ) (0 a T„) 
n n n 

(A.9) =J2J2 C n^k,n^kT n + J2Y,^ zC ^' Bk ^ kTn C nVk,nMdzXn)t T n , 

n fc n fc n k 

where we have used (|A.2|) in the second equality. Now, using (|3.33l) . (|A.8|) and 



{—dt — Afe) Vk, n — —tT n 



dt — A„) 2^ T n — —tT n 



it is easy to show that 



A. 7 Proof of accuracy 

Before establishing our main accuracy result - Theorem 13.41 - we shall need the following lemma. 



Lemma A. 3. Suppose J(y, z) is at most polynomially growing. Then, for every (y, z) and s < t, there exists 
a positive constant C < oo such that for any e < 1 and 5 < I, we have the following inequality 

E y , z [\J{Y s ,Z s )\]<C. 

Proof of Lemma \A . 3[ It is enough to prove the result for J(y, z) = y k and J(y, z) — z k for any fceN. We 
begin with J(y, z) — z k . Under the physical measure P we have 



K[\Z B 



E 



< supE 

5<1 



\zf}\ k ] <C(s,k)<C(t,k), 



by (|A.5|) . Now define an exponential martingale , which relates the dynamics of Z under the risk- neutral 
measure P to its dynamics under the physical measure P. We have 

MP := exp (- jT T(Y S , Z s ) dW s z - \ £ T 2 (Y S ,Z S ) ds^j = ^ . 
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The P-expectation of \Z S \ can be found as follows: 



E 



\Z„ 



= E 
= E 

< |E 



\Z s \ k MP 

\z s \ k exp Q jf r 2 (y„, Z u ) du) (M] 2r ) 



1/2' 



Z a | 2 *expl I T 2 {Y u ,Z u )du 



E 
< ( E 



|z s | 2fe ex P ^ r 2 (r u ,z n )d^ ^ 



1/2 



1/2 



E 



Mi 



2D 



1/2 



(by Cuachy-Schwarz) 
(M( 2r ) is a P-martingale) 



exp( a ||r||^)) <C, 



1/2 



where we have used assumption [3] of section DOl in the last line. We now examine the case J(y, z) — y k . We 
have 



E[|y 5 | fe ] =E 



\Y. 



(i)ife 



< supE 



(i)ife 
A I 



< C(fc), 



by (| A.4|) . Using the same argument as above, one can easily show 



= E 



\2k 



exp (s HAUL)) <C, 



1/2 



)ni*Afj A >] < (Efir^l 2 

which proves the lemma. 

We are now in a position to prove Theorem 13.41 We begin by defining a remainder term R e ' S by 
u e ' 5 = u a 



□ 



o + \/euo,a + VS uo t i + e (u 2 ,o + \/eU3,o) + v 7 ^ (Ve^i,! + e«2,i) + -R e 



The functions uo,Oj ^1,0 and uo,i are the unique solutions to p.30p . (|3.3ip and (|3.32l) respectively. The 
function 112,0 is given by (I3.20[) . And it 2] o is the solution to Poisson equation Q3.13p . To characterize u\ t i 
and tt2,i we must continue the singular perturbation analysis of equation ()3.9j) a bit further. The Q(y/e~5) 
equation that results from continuing the asymptotic analysis is 

(A.10) 0(VeI): = L u 3 ,i + £i«a,i + (-& + £2) + M 3 w 2> o + Mi«i, 



Equation (IA.10[) is a Poisson equation of the form Q3.14p . In order for (jA.lOp to admit a solution 143,1 in 
i 2 (R, 7r), centering condition Q3.15P must in satisfied. In (|A.10p the centering condition corresponds to 



(A.ll) 



= (£i«2,i) + (~d t + (£ 2 » + (M 3 « 2 ,o) + (Mi) u 1>0 . 



Now, by introducing 0) and £(y, z) as solutions to 

^ = / - (/> , 



£ C = r-(r>. 
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and 



and by subtracting (|3.27[) from (|3.26[) . we can express 1*2,1 as 

(A. 12) u 2 ,i = - Qa 2 4>d* x - a-qd^j u ,i - 9 (pxza£,d x - (^d z u ofi + D, 

where D(x, z) is a constant which is independent of y. Substituting (|A.12j) into (|A.11[) characterizes u\ t \ in 
terms of uo,o, M i,0: 1*1,0 and Mo,i- We choose u% t i as the solution to (jA.111) with BC u{0, x, z) = 0. 
Now, we compute 

0= (-dt+L £ ' 5 ) u e - 5 
= {-d t + L € ' s ) R € ' s + -F + -^pF\ +F 2 + VS f-F a + -^F 4 + F 
(A.13) +eRl + VTSR e 2 + 5Rl, 

where 

Fq = LoUo t o, 

F\ = £jQUifl + £iUo,o, 

F 2 = L u 2 ,o + -Citti.o + {-d t + £2) Mo,o, 

F 3 = £ 1io,l; 

F 4 = £ «i,i + £i«o,i + 3Vl 3 u ,o, 

F 5 = £ u 2 ,i + £i«i,i + M 3 ui i0 + Mim ,o + (-<9t + £ 2 ) wo,i, 

i?i = {~d t + £2) u 2fi + £iw 3 , + Ve {-d t + £2) «3,o, 
i? 2 — ( — d t + £2) "i,i + £11*2,1 + Mi«i,o + M 3 m 2 ,o 

+ \/e((-dt +£2)1*2,1 +Mim 2 ,o + M 3 W3,o) +eMiu 3 , , 
i?3 = Miu ,i + M 2 u ,o + M 3 ui,i + (MiUj.,1 + M 2 ui,o + M 3 u 2 ,i) 

+ e(Mitt 2 ,i +M 2 u 2 , ) ■ 

From the choices made in section |3~T1 it is straightforward to show Fq = F\ = F 2 = F 3 = F4 = F5 = 0. 
Hence, from (IA.13[) we have 

(A.14) = {-d t + £ e ' 5 ) i? 6 ' 5 + eR\ + V7S R% + 51%, 

(A.15) R{0,x,y,z) = eG\{x,y,z) + \fe~5 G 2 {x , y , z) , 
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where 

G\(x,y,z) := -u 2 ,o(0, x, y, z) - y^eu 3fi (0,x,y, z), 
G\{x,y,z) := -ui,i(Q, x, y, z) - \feu<i,\(Q, x, y, z). 

Using the Feynman-Kac formula, we can express R e ' S (t,x 7 y,z), which is the solution to PDE (|A.14I) with 
BC (|A. 15|i . as an expectation 

ft 



R e ' (t, x, y, z) = eEa,j/,; 



e -S$Hx.)*»C%(Xt,Y t ,Zt)+ / e" K k(x ^ du Rt,(s, X a , Y s , Z s ) ds 



e-f« k ^ du Rl(s,X s ,Y s ,Z s )ds 



From the assumptions of section IA.3I one can ded uce that the functions (Ri ,R%,Rz,Gt, G%) are b ounded in 



x and at most polynomially growing in (y, z) (see 
by Lemma lA.31 we have 



Fouaue. Papanicolaou. Sircar, and Solna 



mim . Hence, 



\R e > s \ < eC 1 + V7SC 2 + SC 3 < (e + 6) C 4 . 



Finally 



u e ' S - (u ,o + \/e«i,o + uo,i) 

< \R e ' S \ + eu 2 ,o + e 3/ Vo + v / ^ + eVSu 2 ,i 

< (e + S) C 4 + e \u 2 .o + y/eu 3t o\ + VTS \u\ t i + y/eu2,i\ 
<(e + S) C, 



which is the claimed accuracy result. 
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Figure 1: The price of a double-barrier call option is plotted as a function of the underlying x. On the left we 
consider the Black-Scholes model with only a fast- varying factor of volatility Y whose dynamics are given by 
14.31 On the right, we consider the Black-Scholes model with only a slow-varying factor of volatility Z whose 
dynamics a given by 14.41 In each plot, the solid black line corresponds to the full price of the option, the 
dashed line corresponds to our approximation, and the dotted line corresponds to the Black-Scholes price. 
For the plots on the left we use parameters t — 1/12, y = 0, r — 0.05, a — 0.34, p xy — —0.5, f3 = 1, L = 1.5, 
K = 2.0, R = 2.5. For the plots on the right we use parameters t = 1/12, z = 2, r = 0.05, a = 0.34, 
p xz = -0.5, g = 2,L= 1.5, K = 2.0, R = 2.5. 
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Figure 2: The yield of a zero coupon bond is plotted as a function of the time to maturity t. On the left we 
consider the Vasicek model with only a fast-varying factor of volatility Y whose dynamics are given by I4.3I 
On the right, we consider the Vasicek model with only a slow-varying factor of volatility Z whose dynamics 
a given by I4.4I In each plot, the solid black line corresponds to the full yield of the bond, the dashed line 
corresponds to our approximation, and the dotted line corresponds to the Vasicek yield. For the plots on 
the left we use parameters x = 0.03, y = 0, = 0.05, a = 0.02, p xy = -0.5, /3 = 1, O = 0.1 e^ 2 / 4 . For the 
plots on the right we use parameters x = 0.03, z — 1.0, — 0.05, a — 0.02, p xz — —0.5, g = 1, O = 0.1. 
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Figure 3: The implied volatility of a European put option is plotted as a function of the strike price K. 
On the left we consider a JDCEV model with only a fast-varying factor of volatility Y whose dynamics are 
given bv 14.31 On the right, we consider the JDCEV model with only a slow-varying factor of volatility Z 
whose dynamics a given by 14.41 In each plot, the solid black line corresponds to the full implied volatility, 
the dashed line corresponds to our approximation, and the dotted line corresponds to the JDCEV implied 
volatility. For the plots on the left we use parameters t = 1, x = 50, /i = 0.05, a = 10, r] = —1, c = 0.5, 
Pxy — —0.5, y — and (3 — 2. For the plots on the right we use parameters t = 1, x = 50, fi = 0.05, a = 10, 
i] = —1, c = 0.5, p xz — —0.5, z — 2 and g = 2. 
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